Prictional shear cracks 
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We discuss crack propagation along the interface between two dissimilar materials. The crack edge 
separates two states of the interface, "stick" and "slip" . In the slip region we assume that the shear 
stress is proportional to the sliding velocity, i.e. the linear viscous friction law. In this picture 
the static friction appears as the Griffith threshold for crack propagation. We calculate the crack 
velocity as a function of the applied shear stress and find that the main dissipation comes from the 
macroscopic region and is mainly due to the friction at the interface. The relevance of our results 
to recent experiments, Baumberger et al, Phys. Rev. Lett. 88, 075509 (2002), is discussed. 
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A few recent experimental observations of Rubio and 
Galeano 0], and Baumberger, Caroli, and Ronsin Q, on 
the frictional motion of sheared gels sliding along a glass 
surface indicate the existence of self-healing pulses and 
inhomogeneous modes of sliding A regime of peri- 
odic stick slip has been observed in a limited range of 
small shearing rates . It bifurcates towards stationary 
sliding at some critical driving velocity. The slip pulses 
traverse the sample with a velocity much larger than the 
driving velocity but still much smaller than the speed of 
sound. 

Slip pulses in gels seem to be very different from Schal- 
lamach waves and "brittle" pulses studied by Gerde and 
Marder |^ since no observable interface separation oc- 
curs. In this respect, they are more comparable with 
self-healing cracks suggested by Heaton [|| in the con- 
text of seismic events. 

Recent investigations (see, for example, and refer- 
ences therein) point towards an essential importance of 
the underlying friction law in the slip state. It has been 
proved that the simple Coulomb friction leads to the so- 
called "ill-posedness" of the linear stability problem while 
discussing small nonhomogeneous perturbations of the 
stress and strain fields in a sliding mode Moreover, 
Caroli has shown that the existence of slow, periodic 
slip pulses is incompatible with Coulomb friction law. 

In this letter we discuss crack propagation along the 
interface between two dissimilar materials. The crack 
edge separates two states of the interface, "stick" and 
"slip" . We assume that the interface is flat with a strong 
adhesion contact. In principle, we could allow for small 
wavelength surface roughness, but in this case we con- 
sider length scales larger than the longest wavelength 
component. In the presence of roughness, the assump- 
tion of strong adhesion and full contact at the interface 
presumably is only reasonable for "soft" materials with 
relatively small shear modulus. Gels are clearly materials 



of this sort. 

In the slip region we assume a simple linear viscous 
friction law, namely, the shear stress is proportional to 
the sliding velocity. This from the theoretical point of 
view strongly motivated law is usually not discussed in 
literature since it does not lead to the so-called static fric- 
tion phenomenon observed experimentally. However, we 
will see that in our description static friction appears in a 
natural way as the usual Griffith threshold for crack prop- 
agation. The important point is that before the system 
goes into a sliding mode the slip pulse should traverse the 
sample. This requires finite shear stress since the stick 
state of the interface is energetically more favorable. 

With the linear viscous friction law we find conditions 
for crack propagation and calculate the crack velocity as 
a function of an applied shear stress. We find that the 
main dissipation comes from the macroscopic region and 
is due to the friction at the interface. This situation is 
very different from the usual crack propagation where 
the main dissipation is localized in the microscopic tip 
region. 

We also shortly discuss frictional shear cracks inside 
homogeneous materials. The point here is that in mode II 
(and in mode III) cracks there is no macroscopic opening. 
If two surfaces remain in contact, the standard bound- 
ary conditions, namely the vanishing of normal and shear 
stresses on the crack surfaces, are not theoretically mo- 
tivated. The relative sliding velocity of the two surfaces 
should lead to nonzero shear stresses. Finally we discuss 
the relevance of our results to the experimental observa- 
tions II . 

Consider an elastic solid sliding on a flat rigid sub- 
strate. Assume that the elastic solid occupies the space 
H > y > 0, and let (x,y,z) be a coordinate system 
with the plane y — corresponding to the surface of the 
solid, see Fig.l. We discuss the plane strain situation 
with tt^ = 0, where u is the displacement vector. We 
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assume that the interface can be in two states: "stick" 
and "sHp". The boundary between these two states is 
described by the crack edge which moves with a veloc- 
ity Vtip in the x direction. In the stick region the dis- 
placements are continuous and since we assume a rigid 
substrate, the boundary conditions are: Ux = Uy — for 
X — Vtipt > and y = 0. In the slip region we assume 
that the two solids (for all times) are in contact, Uy — 
for X — Vtipt < and y = 0, while we allow for a finite 
relative sliding velocity iix- This sliding velocity leads to 
frictional shear stress at the interface where we assume a 
linear viscous friction law: 

axy = aiix, (1) 
with a being the viscous friction coefficient. 




FIG. 1. An elastic body sliding on a rigid substrate 



It is reasonable to assume that the interface energy in 
the stick phase is smaller than the interface energy in the 
slip phase since the adhesion contact in the stick region is 
stronger. Let us denote this energy difference by 7. It is 
clear that without external loading the stick phase is en- 
ergetically favorable and a finite shear stress is required 
to get the interface into the slip state. Let us assume that 
far ahead of the crack tip, the solid is homogeneously 
strained with and stressed with Uxy = where 
/i is the shear modulus. The strain energy is ^^{u'^yYH. 
Far behind the crack tip where the stress is relaxed, only 
the interface energy 7 remains. The slip state will be 
realized only if A = /i(u!^)^iJ/7 > 1. In this case the 
crack should propagate in the positive x direction. Oth- 
erwise, the crack would propagate with negative velocity 
and the stick phase will be restored. Condition A = 1 is 
nothing but the usual Griffith threshold for crack propa- 
gation. On the other hand, in the context of the friction 
problem, this condition may be interpreted as a static 
friction threshold: a finite shear loading is required to 
get the system into the sliding mode. 

If the whole interface is in the slip state, a steady-state 
motion of the elastic body is possible with a velocity 

V = 2Suxy, (2) 

where S = fi/a is the velocity scale given by the fric- 
tion law. We note that this homogeneous sliding mode 



is linearly stable for any velocity with respect to small 
nonhomogeneous perturbations of the stress and strain 
fields localized in the surface region. In this respect the 
viscous friction law is very different from Coulomb fric- 
tion which leads to a linear instability and ill-posedness 
of the problem as it has been intensively discussed in 
literature 

On the other hand, the homogeneous sliding mode may 
be unstable against the resticking pulse (nonlinear "heal- 
ing instability") if the corresponding value of A < 1. 
Since in this case, the strain which defines the value of A 
is related to the steady-state sliding velocity by Eq.(2), 
we find that the homogeneous sliding is stable against the 
healing instability only above the critical sliding velocity 

K = 25(7//ii/)'/'. (3) 

Now let us turn to the calculation of the crack tip ve- 
locity Vtip as a function of the dimensionless driving force 
A. The strategy is as follows: we solve the elastic prob- 
lem in the vicinity of the crack tip and then calculate 
the energy fiux into the crack tip and the dissipation due 
to the friction at the interface. Finally, using the energy 
balance, we find the crack velocity. 

Let us start from some qualitative estimates. Assume 
that, as in the usual crack problem, the singular behav- 
ior of the displacement vector is given by a square-root 
singularity in the vicinity of the crack tip. Then the dis- 
sipation rate at the interface 

Jd^a J {iixfdx (4) 

diverges logarithmically. This is already quite a remark- 
able observation. In the usual crack problem the main 
dissipation comes from the close vicinity of the tip and of- 
ten requires the introduction of microscopic models. Here 
we have the chance to avoid such a detailed microscopic 
description by using some microscopic length scale as a 
cutoff which enters only the logarithm in the final result. 
Note that dissipation due to the bulk viscoclasticity ef- 
fectively only leads to tip dissipation on microscales. In- 
deed, the viscoclasticity gives corrections to the stress 
tensor of the form "quik- The dissipation rate due to this 
effect diverges strongly at small distances (as 1/r) and 
correspondingly decays at macroscopic distances. Thus, 
this effect can be incorporated into the tip dissipation. 

Now we solve the elastic problem more accurately, 
while still using a quasistatic approximation for the mo- 
ment. The generalization to the full elastodynamic de- 
scription is straightforward and will be given below. In 
the co-moving frame of reference and in the vicinity of the 
crack tip, the displacement field for the static elasticity 
and the boundary conditions formulated above is 

Ux = ARe [yz^^^ - i (3 - 4i/) z^/X] , 
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ARe [ly. 



(5) 



Here z is a complex coordinate z — x + iy, h' is the Pois- 
son ratio, and ^ is a real amplitude. The spectrum of A 
is purely real and given by the following equation: 



exp (i27rA) = — 



l + i7re/2 
1 - iTTe/2 ' 



with 
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In the limit of small values of |e| for the leading crack 
displacement component, we have 2A ~ 1 + e. 

Having defined the displacement field we can calculate 
the energy flux = cTikUk and the local energy release 
into a small semicircular region with some microscopic 
radius a around the crack tip, 
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Jo = 27r^(3 - Av){l - v)VupA^a: 



(8) 



The dissipation rate due to the interface friction with ex- 
clusion of the small region of size a close to the tip is 
given by Eq. (4): 



Jrf = a(3 - AvfA^V^^^e-^ {h^ - a' ] = Jq 



{H/af 1 
(9) 



Here H = with / being an undetermined yet 

function of the Poisson ratio; _ff is a thickness of the 
sample. This function can be found by solving the elastic 
problem for given geometry and all boundary conditions. 
As we will see the tip velocity does not crucially depends 
on the actual value of the factor /, which is of the order 
of unity. 

On the other hand the local energy release into the 
crack tip Jo must compensate the surface energy differ- 
ence: Jo = ^Vtip. Note, that here we have neglected 
the dissipation at the tip in comparison with the energy 
release "fVup, which is reasonable at small tip velocities 
compared to the velocity of sound. Finally, using the 
global energy conservation law. 



Joo = Jo + Jd = M (u^j) HVtip, 



we find 
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Since e is given by Eq. (7) , this result is a compact repre- 
sentation of the crack velocity as a function of the driving 
force A — n (w^y) H/j. Note that A = 1 corresponds 
to Griffith equilibrium. Eq. (11) is valid for small |e|. 
The explicit expression for the crack velocity reads 
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In the case of small A — 1 we obtain 
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This result corresponds to a small dissipation rate com- 
pared to the total energy fiux, Jd <^ Joo , and can also be 
obtained using perturbation theory: we solve the elastic 
problem neglecting friction (axy = at the interface) and 
then calculate the dissipation rate (4) using this solution 
as zero order displacement field. 

Up to now we have used the static approximation. The 
tip velocity should be small compared to the sound veloc- 
ity. The elastodynamic generalization is straightforward. 
Using the standard approach to the singular solutions 
near the crack tip H , we find the displacement field u 
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with = 1 



{Vup/Cdf and al = \- {Vup/C,)\ where 
Cd and Cg are dilation and shear wave speed. The spec- 
trum of A is still given by Eq. (6) but now e reads 
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Thus, elastodynamic effects just lead to a redefinition of 
£ in the main result, Eq. (11), which remains valid. For 
small velocities, Eqs. (14) and (15) reduce to Eqs.(5) and 
(7), respectively. 

The most serious problem with large velocities arises 
in connection with a self-consistent description of the dis- 
sipation at the tip. This part of kinetics cannot be con- 
sidered macroscopically for arbitrary tip velocities. One 
can only treat the case of small velocities, Vtip ^ Cs in 
a model independent way, by introducing the tip kinetic 
coefficient. For higher velocities the so-called velocity de- 
pendent fracture energy ^{Vnp) is introduced. This func- 
tion contains information about the usual surface energy 
7 and tip dissipation and reduces to the surface energy 
in the static limit. The main dissipation in our approach 
arises from the friction between both sides of the crack 
and can be treated macroscopically. Note that this part 
of the dissipation can be described by the velocity inde- 
pendent friction coefficient a even at large tip velocities 
in the case of small sliding velocities. 

Up to now we have discussed shear cracks along the 
interface between two dissimilar materials (the case of a 
rigid substrate). We also present the result for the case of 
frictional shear cracks inside homogeneous elastic materi- 
als. Such cracks can propagate under the shear loading in 
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amorphous materials, or along grain boundaries in crys- 
tals (in single crystals the well known dislocation mech- 
anism of plasticity should be favorable). The boundary 
conditions on the crack surfaces, which remain in con- 
tact, are: continuous normal displacement and continu- 
ous normal and shear stresses. The shear stress is also 
given by Eq.(l) with iix being the relative sliding velocity 
of two crack surfaces. We note that these boundary con- 
ditions are quite different from the standard boundary 
conditions of mode II (and III) cracks: zero normal and 
shear stresses on the crack surfaces In our case, the 
frictional shear crack, we find the following expression for 
e which enters the general result, Eq.(ll): 

««(!-«') Vup . . 

Now let us discuss the relation of our results to exper- 
imental observations of Baumberger, Caroli and Ronsin 
They performed experiments of a gel sliding on a 
glass plate. The driving velocity was given and the shear 
force and thus the average shear stress was deduced from 
the spring elongation. Above some critical driving ve- 
locity Vc ~ 125fj.m/s, steady sliding was observed. At 
velocities smaller than the critical one, periodic stick slip 
sets up (see figures in Q). Upon increasing the driving 
velocity V, no hysteresis of the transition was detected. 
In the stick slip regime they observed the propagation 
of self-healing pulses with no opening, nucleated periodi- 
cally at the trailing edge of the sample. The propagation 
velocity of these cracks was about 60 times larger than 
the critical sliding velocity, yet still much smaller than 
the shear wave speed. 

The existence of a critical sliding velocity, where sta- 
tionary sliding is stable against the healing instability, 
appears naturally in our description and is given by Eq. 
(3). The characteristic value of the shear strain in the 
sliding mode near the critical velocity experimentally 
was about u^y = 0.04. Thus, we can estimate from 
Eq. (2) the characteristic velocity S = l.bmm/ s and 
from Eq.(3) we find that the characteristic difference be- 
tween the interface energies in the slip and stick states is 
7 = O.IJ /m?. One would expect that for ordinary elas- 
tic materials the velocity 5" should be of the order of the 
speed of sound. However, for gels the shear modulus /i 
is much smaller than for ordinary materials. The shear 
wave speed Cg = (/i/p)^/^ is only 2m/ s. The velocity 
S — /i/a is linear in /i and should be even smaller. This 
is a possible explanation for a relatively small value of S 
compared to Cs- 

In the stick-slip regime, which exists below 14, the nu- 
cleation of a slip pulse takes place at the trailing edge of 
the sample and requires overshooting above the Griffith 
threshold according to the experiment. This overshoot- 
ing is not small and in order to estimate the crack veloc- 
ity we can use Eq. (12) since the velocity is still much 



smaller than the speed of sound. Because of the week log- 
arithmic parameter dependence, we conclude that Vup is 
of the order of S and essentially independent of the driv- 
ing velocity in agreement with experimental observations. 
After the slip pulse traversed the sample the stress drops 
below the Griffith threshold and resticking takes place 
via propagation of a healing pulse. Its velocity is still de- 
scribed by Eq. (12) with A < 1. This periodic stick slip 
regime bifurcates towards stationary sliding dX V — Vc 
where stresses are always above the Griffith threshold. 
For driving velocities slightly below Vc , characteristic val- 
ues of A for resticking are close to 1. Since the velocity 
of the resticking crack is small and comparable with the 
driving velocity in this range, a complicated collective 
behavior of self-healing pulses is observed |^ . 

While our results are in a qualitative agreement with 
experimental observations, we still underestimate the 
crack velocity which is a few times smaller than in the 
experiment. The other discrepancy is due to the ob- 
served nonlinear behavior of the stress with the veloc- 
ity in the stationary sliding regime (the so-called shear- 
thinning rheology). On the other hand, the geometry of 
the experiment is such that the total macroscopic friction 
of the sliding sample obviously depends on the processes 
taking place at the edges of the sample. The stresses 
here are highly inhomogeneous and the kinetic phenom- 
ena should be considered with a great care. We think 
that these boundary effects may be responsible for the 
discrepancies mentioned above. 

Another subject of future investigations should be the 
collective behavior of self-healing pulses in the spirit of 
Ref. [Q. Further theoretical and experimental investiga- 
tions are needed to shed light on this phenomenon where 
two intriguing problems, crack propagation and friction, 
combine together. 

Discussions with S.V. lordanskii, B.N.J. Pcrsson and 
D.E. Temkin are greatly appreciated. V.M. thanks 
Forschungszentrum Jiilich for its hospitality. 
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